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A discrete-time Quantum Walk (QW) is essentially a unitary operator driving the evolution of a single particle 
on the lattice. Some QWs admit a continuum limit, leading to familiar PDFs (e.g. the Dirac equation). In this 
paper, we study the continuum limit of a wide class of QWs, and show that it leads to an entire class of PDFs, 
encompassing the Hamiltonian form of the massive Dirac equation in (1 + 1) curved spacetime. Therefore a 
certain QW, which we make explicit, provides us with a unitary discrete toy model of a test particle in curved 
spacetime, in spite of the fixed background lattice. Mathematically we have introduced two novel ingredients for 
taking the continuum limit of a QW, but which apply to any quantum cellular automata: encoding and grouping. 
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I. INTRODUCTION 


Quantum walks (QW) were originally introduced lUSl as 
dynamics having the following features: (i) the underlying 
spacetime is a discrete grid; (ii) the evolution is unitary; 
(in) it is homogeneous, i.e. translation-invariant and time- 
independent, and (iv) it is causal, i.e. information propagates 
strictly at a bounded speed. 

Quantum Computing has a number of algorithms that are 
phrased in terms of QWs, see 0 for a review. Our focus here 
is on QWs models per se, or as models of a given quantum 
physical phenomena, through a continuum limit. Such QWs 
models have a broad scope of applications: 

• they provide quantum algorithms, for the efficient simu¬ 
lation of the modelled phenomena upon a quantum sim¬ 
ulation device 0; 

• even for a classical computer they provide a stable nu¬ 
merical scheme, thereby guaranteeing convergence of 
the simulation as soon as the scheme is consistent 0; 

• they provide discrete toy models to explore founda¬ 
tional questions iHni. 

In this paper we introduce Paired QWs, which are both a 
subclass of the general QWs described above, and generaliza¬ 
tion of the most usual QWs found in the literature. Basically, 
(i) the input is allowed a simple prior encoding and (ii) the lo¬ 
cal unitary ‘coin’ is allowed to act on larger than usual neigh¬ 
bourhoods. Moreover, the coin is allowed to depend on space 
and time, as in other QW models. 

We show that Paired QWs admit as continuum limit the 
class of PDFs of form 

dtipit, x) = x) + x) + iCipit, x) (1) 



FIG. 1: Usual QWs. Times goes upwards. Fach site contains a 2d- 
dimensional vector 'ip = tp^ © ip~. Fach wire propagates the d- 
dimensional vector ip^. These interact via the 2d x 2d unitary W. 
The circuit repeats infinitely across space and time. Notice that there 
are two light-like lattices evolving independently. 


with Bi and C hermitian and |i?i| < Id. This class of PDFs 
includes the Hamiltonian form of the massive curved Dirac 
equation in (1 + 1)-dimensions lfl4l for any bounded metric 
in any coordinate system, together with an electromagnetic 
field. Given the PDF we wish to simulate, we are able to 
retro-engineer the corresponding Paired QW. 

The results extend the connection between QWs and the 
Dirac equation, first explored in iHHaiia, and further devel¬ 
oped in EllElIIMIll. Fxtension to curved spacetime was 
initiated in ll20l - [22l . more carefully discussed in the conclu¬ 
sion. 

We proceed by first formally defining the model, in Section 
[n| We then compute the conditions for the continuum limit to 
exist, and provide a complete parametrization of the QW op¬ 
erators in terms of the metric, in Sectionjl^ Then we identify 
the continuum limit with the Dirac equation in curved space- 
time, and validate the model with numerical simulations in 
Section [TV] F inally, we discuss perspectives and related works 
in SectioniVI 
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II. MODEL DEFINITION 


Usual one-dimensional QWs act on the space © 

C'^), equal to the set of square summable sequences in the 
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FIG. 2: The input to a Paired QW is allowed to be encoded via a 
unitary E, and eventually decoded with E"' . 



FIG. 3: When the scheme is iterated, the decoding of the previous 
time-step cancels out with the encoding of the next time step. Thus 
the only relevant encoding/decoding are those of the initial input and 
final output. A Paired QW is therefore really just a QW, with a par¬ 
ticular choice of initial conditions. 

space 0 C*^). Often, the dimension of the internal 

degree of freedom is two, corresponding to d = 1. We de¬ 
note those functions taking a lattice position x into the 
^2d-vector (t, x)(Bip~ {t,x), where each (t, x) is a 
vector. 

These QWs are induced by a local unitary W from 
to often referred to as the coin. Hence c = 2d k 
often referred to as the coin dimension or internal degree 
of freedom of the walker. The reason why c must split as 
d + dis because of the way W is wired: each W{t, x) takes 
one half of ip{t,x — 1) (more precisely, its d upper compo¬ 
nents X — 1)) and half of tpit, a; 0 1) (more precisely, 
its d lower components ip~{t,x + 1)) in order to produce 
+ l,a:). This way the inputs and outputs of the differ¬ 
ent W{t, x) are non-overlapping and they can be applied syn¬ 
chronously to generate the QW evolution over the full line, so 
that 

U{t)-.= ^W{t,x) (2) 

generates one time step of the QW (we remark that t indicates 
possible time dependence of the local unitaries; do not con¬ 
fuse U{t) with the evolution operator from time 0 to time t). 

It follows that usual QWs evolve two independent light¬ 
like lattices, as emphasized in Fig. On one of the light-like 
lattices, the evolution is given by 

Vit) := 0 W{t, x) and V{t + 1) := 0 W{t + l,x). 

xe2Z xe2Z+i 

(3) 

whilst on the other lattice everything is shifted by 1 in posi- 



FIG. 4: A Paired QW obtained by spacetime grouping of an ordinary 
QW. The green triangles define the appropriate encoding E, that re¬ 
lates the fine-grained input ip with the coarsegrained input tp'. The 
dotted line indicates at + 2 space-like surface fine-grained output. 
This surface is recovered by undoing the triangles above these dotted 
line, which is the role of EK 


tion, 

V{t):= 0 W{t, x) and V{t + 1) := 0 W(t + l,x). 
xe2’z+i xe2Z 

(4) 

Paired QWs arise as follows. Bunching up every 'ippt, x—1) 
and tp{t, a: 0 1) site into (p{t, x) = tp{t, x — 1) 0 tp{t, x 0 1), 
and applying a unitary encoding E to each bunch, we obtain 
(p'{t, x) = E(p{t, x). We may now define a QW over the space 
these encoded bunches cp'. The local uni¬ 
tary W' will be from to and each W'{t, x) will take 
one half of (p'{t, x — 2) (more precisely, its 2d upper compo¬ 
nents) and half of (p{t, x + 2) (more precisely, its 2d lower 
components) in order to produce (p'(t + 2,x). The inputs and 
outputs of the different W'{t,x) are again non-overlapping 
and they can be applied synchronously to generate the QW 
evolution over the full line, 

U{t):= ^W\t,x). (5) 

In the end, each (p'{t+2, x) maybe decoded as (p{t+2, x) = 
E^'{t 0 2, x) and be reinterpreted as (p{t 0 2, x) = ip{t 0 
2,x—l)0V'(f02,x0l). Clearly this Paired QW (pictured 
in Figs. [^and[^ phrased in terms of cp' and d! = 2d is no 
different from the usual QW definition right above. At least 
from a discrete point of view. 

When looking for a continuum limit, a subtle difference 
arises. Indeed, say that the regular initial condition is given 
in terms of the fine-grained spacelike surface of tpit), which 
is assumed to be smooth, i.e. ipit, x) ~ tp(t, x + 1). Then the 
resulting (p{t) will be smooth both externally, i.e. (p{t, x) « 
(/)(<, X 0 1), and internally, i.e. (p{t,x) « ip{t,x) 0 ip{t,x), 
which is not so usual to ask for. Similarly, (p'(t) will be 
smooth both externally, i.e. (p'{t,x) « (p'{t,x 0 1) and in¬ 
ternally, <p'{t,x) « E{'ip{t,x) 0 ip{t,x)). It turns out that 
such reinforced regularity conditions are necessary for some 
Paired QWs to have a limit. 

The next paragraph is to emphasize that Paired QWs, and 
their reinforced regularity assumptions, are not ad-hoc: they 
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arise naturally when one performs spacetime grouping of 
QWs. 

A natural example of Paired QW is provided by performing 
the spacetime grouping of a usual QW, an operation which we 
now explain. The spacetime grouping operation takes a QW 
over 0z(C‘^ 0 C^), with local unitary W into a Paired QW 
over with local unitary W, as pictured in 

Fig-0 

It is important to notice that if the initial condition was 
given by for the original walk, the initial condition for 
the spacetime grouped QW is now given by the = 

E{x)(j){t, x), and x) = a; — 1) © x + 1), as pic¬ 
tured in Fig. 0 again. In the end, each + 2,x) may be 
decoded as (pit + 2,x) = {x)p'{t + 2, x) and be reinter¬ 

preted as p{t + 2,x) = tp{t + 2, a: — 1) © ^{t + 2, a; + 1). 
This spacetime grouping is reminiscent of the “stroboscopic” 
approach of Col ED, but has the advantage of mapping usual 
QWs into usual QWs of increased dimension. 

This paper studies the continuum limits of Paired QWs (not 
necessarily arising from a spacetime grouping) for d — 1, 
systematically. Recall that for d = 1: ip{t) is in £^(Z;C^) 
and represents the ‘physical’ field; p'{t) is in f^(Z;C‘*) and 
represents a paired, encoded version of it; W' is the 4x4 
coin operator. For our purpose, it will be useful to redefine the 
bunching-up p{t,x) as 


with 


(p{t,x) 


'u{t, x )' 
d{t, x) 
u'{t,x) ’ 




d'it, x)_ 

u(f, x) 

= H 

p~^it,x + 1) 

u’it, x) 

_p+it,X- 1)_ 

dit, x) 

= H 

'p~it,X + 1) 

d'it, x) 

tp-it,x- 1)_ 


( 6 ) 


(7a) 

(7b) 


where El 


1 

©2 


1 

1 



is the Hadamard matrix. 


Notice that the ip{t, x) dependencies are the same as before, 
this is just a matter of applying a unitary pre-encoding. This 
convenient choice of basis is so that in the continuum limit, 
u{t,x) becomes proportional to ip~^{t,x), whereas u'{t,x) 
becomes proportional to the spatial derivative of ip^{t, x). 

Armed with those conventions on Paired QW, we can focus 
on how pout '■= -\- 2, x) gets computed, from pin := 

pit, X — 2.)(Bpit,x + 2). This ©C^ to function may be 
thought of as the local rule of a cellular automata with cells in 
Its explicit formula is given by 


G = E\t + 2, x)W'it, x)iP' © P) 

{Eit,x - 2) ® Eit,x+ 2)), ( 8 ) 


where the 2x4 projectors P and P' pick-up the u, d (non- 
primed subspace) and u',d' (primed subspace) coordinates, 
respectively. Thus 


p{t + 2,x) = Gip{t, X — 2) ® Pit, X + 2)). (9) 


III. CONTINUUM LIMIT 


From now on, we consider that t and x are continuous vari¬ 
ables, and choose the same discretization step e G M’*' for 
each coordinate. In particular note that, to first order in e, 
we have that u ~ y/2'p'^, d ~ '/2p~, u' ~ e^/2dx'p'^ and 
d' ~ £y/2dxp~ see 0 - To start investigating the contin¬ 
uum limit of the system defined by Eq. 0, we compute the 
expressions for the input and output. 

The expansion of the input to first order in e in terms of 
u, u', d, d' is 


'u 


'u 


■-2m'- 


■2m'- 

d 


d 


-2d' 


2d' 

0 

© 

0 

+ 

u' 

© 

u' 

.0. 


.0. 


d' 


d' 


( 10 ) 


We stress that u' and d' are themselves proportional to e, 
hence the last term is proportional to e. 

The expansion of the output to first order in e in terms of 

u, u', d, d' is 



'u 


'2£dtu 

4^out (j^: — 

d 

0 

+ 

2£dtd 

u' 


.0. 


d' 


( 11 ) 


Next we specify the structure of the walk and encoding op¬ 
erators. We shall assume, for simplicity, that the matrix ele¬ 
ments of W and E are analytic functions of (f, x) and e. 

First, we set W' := with unitary and W 

hermitian. This enforces the unitarity of W, and is without 
loss of generality, since only its expansion to first order in e 
matters: 

Wit,x) ~ W^^\t,x) + i£W^^\t,x)Wit,x). (12) 

Then, in a similar manner, we define E := with 

unitary and E hermitian. Hence, to first order in e, 

Eit,x) ~ E^^\t,x) +i£E^°\t,x)Eit,x). (13) 

Here is some notation. Any matrix A G will be writ¬ 
ten in block form as A = ^ 

j = 1,..., 4. Let X = Ux ® I,Y = ay ® I and Z = Oz® I, 
where (ct^, Uy, az) are the Pauli spin matrices. 

Notice that, for any A G the following simplifica¬ 

tions hold: 

iP'® P)iA® A)iv ® v) = XAv Vu G (14) 


and 

iP'® P)iA® A)i-v®v) = XZAv VuGC^. (15) 

Next we develop the zeroth order and the first order expan¬ 
sion in £ of Eq. 0. 
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A. Zeroth order 

For the left hand side we have just the zeroth order of o, 
while for the right hand side there is only one term which 
does not contain e, obtained multiplying all the zeroth order 
contributions. Hence 


~u 

d 

0 

= F;(o)tpp(o)(p/ 0 p)y(o) © i;(o)) 

'u 

d 

0 

© 

'u 

d 

0 

.0. 


.0. 


.0. 




u 

d 

0 ’ 
. 0 . 


(16) 


with 


B = E^°'>'<ZE^°'> 

(21a) 

N = {dtE^°^^)E^°'' 

(21b) 


(21c) 


(21d) 

} we shall study separately what happens in 


the primed and in the non-primed subspaces. 


C. Continuum limit equation 

Projecting Eq. ( |20| ) on the non-primed subspace, we obtain 
an equation with time derivatives in the left hand side, 


where we used the simplification ( [l4l i. The only non-trivial 
relations are 


2edtu 

2edtd 


= B, 


2u'' 

2d' 


+ e(2iVi + Tj -f 2Mi) 



( 22 ) 


u 

d 

o' 

0 


^^(0)tpp(0)jf^(0)^ 


u 

d 


^^(0)tpp(0)jf^(0)^ 


u 

d 


(17) 

(18) 


To satisfy ( [T7| ) for arbitrary u and d, we must take the identity 
for block 1. Now, since the matrix in ( |T6] l is unitary, then both 
its rows and its columns must sum to one, thus the blocks 2 
and 3 become zero, and ( fTS) ! is automatically satisfied; we are 
left with the choice of an arbitrary unitary U G U{2) for block 
4, to complete the matrix. Hence 


i;(o)tvp(o)xF;(o) = J © [/, (19) 


where the direct sum is respect to the non-primed subspace 
(spanned by the first two entries) and the primed subspace 
(spanned by the last two entries). 


B. First order 

For the left hand side we have just the first order of 
note it contains time and space derivatives of tjj^. For the right 
hand side we multiply and collect all possible combinations 
in which only one term contains e. Then, after a long but 
straightforward calculation (see Appendix [A] for the details), 
we get 


'2edtu 


'O' 


-2u'- 

2edtd 

u' 

= (/©(/) 

0 

u' 

+ {I(BU)B 

2d' 

0 

d' 


d' 


. 0 . 


+ e\^i2N-iE){I®U) 
+{I®U){iE + 2M) + T^ 


u 

d 

0 

. 0 . 


( 20 ) 


Switching to the original coordinates, and writing 

= [^+it,x),ip-{t,x)]^, 

dttpit, x) = x) + ^A^i + y + x). 

(23) 

From pia|i, applying Leibniz rule and using (|21d|i we 


have ll29ll 

d^B = M + = 2^M. (24) 

From ( |21b| ), the unitarity of implies that N is skew- 
hermitian, 

= -N. (25) 

From ( |21c| i, 

T = (26) 

= ( 21 ) 

= i(7 © U)E^°'>'^XWXE^°\ (28) 


where we used the zeroth order condition ( |l9] l. Inverting, 

i£;(o)txiEXF;(°) = {I(BU^)T= . (29) 

Since the left hand term is skew-hermitian we have that 

tI = -Ti (30a) 

T3 = -TIu (30b) 

tIu = -U^T4. ( 30 c ) 

Therefore, by spliting Mi into its hermitian and skew- 
hermitian parts, and using equations ( |24l i, ( |25] l and ( |30a| i, the 
continuum limit has the general form 

dtip{t,x) = Bidxi’{t,x) + ^dxBiip(t,x)+iCilj{t,x). (31) 

where C is an hermitian matrix defined by 

iC = A^I + y + i3Ml. 


(32) 
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D. Compatibility constraints 


2. Determination of and 


Projecting Eq. onto the primed subspace, gives 


= U 


2 UB 2 


+ £ ( 2N2 — '1E2 


+ iUE 2 + 2UM2 + T2 


(33) 


Eq. ( [33l l does not involve time derivatives. Therefore, these 
equations must be understood as constraints. Where do these 
come from? Recall that the aimed continuum limit equation 
Q is over a field, but the QW employed is over the field 
obtained by pairing it. Thus, the C'* field has some internal 
smoothness initially, which the QW must preserve. More pre¬ 
cisely, in order to have nontrivial, time-dependent solutions, 
the coeffients of [it, and [it', v'Y must vanish separately: 

r U{IE2B2) = I, (34a) 

1 2 N 2 - i(/ - U)E2 + 2 UM 2 + T2 = 0. (34b) 


Erom Eq. ( |2Tal i we know that diagonalizes B. Then, 
its columns can be chosen to be any complete set of nor¬ 
malized eigenvectors of B. More generally, we could take 

for arbitrary R,Sg U{2), because of the de¬ 
generacy of order two for each eigenvalue -fl, —1. 

Eor the special case of B in Eq. ( [36l l, we can give an explicit 
solution e‘'°\ 


E^^^ = 



(4 

0 


0 \ 

1 

0 


0 

1 

to 

71 


0 


0 



^2 

0 

to 


(37) 


where vf" = y/1 ± di, i G {1,2}. 

Once E^^'> is known, we can compute from ([T9]l. 


3. Determination of E and W 


E. Existence of solutions 

Up to now we have determined the continuum limit, pro¬ 
vided that the constraints ( |34al i-( |Mb] i are satisfied. In this sec¬ 
tion we show that, given any hermitian Bi and C, there are 
indeed compatible choices of W and E. 

The strategy is the following. Eirst we show that Bi along 
with constraint ( |34a| l determines the zeroth order part of E 
and W'. Then, using C and ( |34b| l we complete the solution 
by determining the first order terms. 


The choice of E^^'> determines Ni and Mi via Eqs. ( |2Tbl ) 
and ( |21d| l, so Ti is fixed by Eq. ( [3^ once we choose C. 

Since E does not appear in the continuum limit, without 
loss of generality we can take i? = 0. In this way T 2 is hxed 
by the contraint ( |34b| i. 

In order to complete T it is now sufficient to take = 0, 
and T3 from ( |30b| i. Einally, from ( |29l l we hnd W, 

W = -iXE^°^ (J 0 UYtE^°^X. (38) 


F. Recap 


1. Determination of B and U 

Consider the spectral decomposition Bi = VDV\ D = 
diagjdi, d2}- In Appendix | b] we show that Eq. ( |21al i implies 
that di, d 2 must belong to the interval [—1,1], and provide the 
general form of B given Bi (see section [TV] for a discussion 
about the eigenvalue constraint). Here we just pick one par¬ 
ticular solution, namely 

® = (‘'o' F.) B (0 ?) • 05) 

where B is 


We have shown that the continuum limit of our model is 
given by Eq. Q. Moreover, we provided a procedure to ob¬ 
tain the parameters of the quantum walk, namely the unitaries 
W' and E, given a pair of hermitian matrices Bi and C, pos¬ 
sibly spacetime dependent. We remark that the choices made 
in the procedure are in general not unique. In the emergent 
continuum limit, different choices of E, W' lead in general to 
the same equation. 

Notice also that the minimal coupling (e.g. electric field) is 
already considered in the parameter C. 

The whole procedure was programmed in sagemath, and 
made available in li2^ . 


B = 


di 

0 

-Aie-*''! 

0 


0 

d2 

0 

-A2e-‘''^ 


-Aie*''! 

0 

-di 

0 


0 

-A2e*''= 

0 

—d2 


(36) 


with A; = Yl — di, sin Pi = ±|cii|, — 7r/2 < rji < 7r/2, 
i G {1,2}. 

Note that U is now fixed by Eq. (|34a|). 


IV. RECOVERING THE DIRAC EQUATION 


On a spacetime with metric tensor and in the absence 
of external helds, the Dirac equation in Hamiltonian form iflTll 
is idttp = Huf), with 


Hd = -i 



dx - ^dx 



-P 
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FIG. 5: Simulation of the Paired QW for the Schwarzschild met¬ 
ric with mass parameter M = 0.5. We plot the probability density 
for a particle with initial condition given by a gaussian wavepacket 
4>{x) oc / (r6+(p)-F M_ (p)) dp where 

a;o = 3.0, po = 50, o = 1.56 and u± are the eigenvectors of 
the free Dirac Hamiltonian Hq = ap mj3. The mass of the par¬ 
ticle is m = 50. For comparison, we show in yellow a grid of null 
geodesics. The lattice spacing is £ = 5 x 10“®. 


where m is the mass and a,/? are matrices satisfying = 
0^ = I, a/3 + /3a = 0 (here and in the following we assume 
natural units, h = c = G = 1). Here e'^{t,x) are the dyads, 
which are related to the metric via x)e^{t, x)e'^{t, x) = 
r]ab, where r]ab is the Minkowski metric. 

Making the identification with the continuum limit of our 
discrete model, Eq. ([T]i, we find 

Si=-4a-eo (39) 

Sq 

777 

C = -3/3. (40) 

These equations allow to find the QW parameters, associated 
to a given metric. The constraint that the eigenvalues of Bi 
are di^2 G [—1,1] represents the finite speed of propagation 
on the lattice. In practice, for any region of spacetime where 
the metric field is bounded, it is possible to rescale the coordi¬ 
nates in such a way that the physical lightcones are inside the 
“causal lightcones” of the discrete model. 

For instance, we can specialize the previous considerations 
to the case of the Schwarzschild metric, whose radial part is 

ds^ = (1 — 2M/x)d0 — (1 — 2M/x)~^dx^, (41) 

where x corresponds to the radial coordinate l30l- Choosing 
the chiral representation namely a = cr^ and P = ax, 


we have 


Bi 






c = 




(42) 

(43) 


A simulation for a particle in the Schwarzschild metric is 
shown in Fig. Again the sagemath program which con¬ 
verts any metric into the corresponding QW, and produces 
such simulations, is available at IJS). 


V. SUMMARY, RELATED AND FUTURE WORKS 

In summary, we have constructed a QW, i.e. a strictly uni¬ 
tary, causal, local evolution, which implements the idea of 
“encoding-evolution-decoding” of the discrete dynamics, de¬ 
fined in Eq. We have found that Paired QWs with d = 1 
admit as continuum limit all PDFs of the form Q. This class 
encompasses the Hamiltonian form of the Dirac equation in 
curved spacetime, together with an electromagnetic field. 

In this model, curvature is effectively implemented by a 
collection of spacetime dependent local unitaries (namely, the 
unitary operators E{t, x) and Wit, x)), which are distributed 
over a fixed background lattice, and whose purpose is to drive 
the particle according to the metric. 

Our main contribution in this paper is to extend the existent 
literature in two directions: we allow for massive particles and 
for any bounded metric field over an arbitrary coordinate sys¬ 
tem. Let us explain the relation between the present paper and 
related works II20II2TI in more detail. In Il20ll2n . Di Molfetta 
et al. have systematically studied the continuum limit of a 
two-time-step QW with arbitrary coin. This stroboscopic ap¬ 
proach made it possible to recover the (1 -F 1)-Weyl equation 
in curved spacetime, for metrics of having poo = 1- Such 
metrics exclude the Schwarzschild metric of Fig. for in¬ 
stance. Yet, any metric can be brought to have poo = 1 if one 
is allowed a prior change of coordinates. On the other hand, 
in (1 -F 1)-dimensions, every space-time is conformally flat 
up to a change of coordinates, and therefore allowing changes 
of coordinates could be said to oversimplify the problem (for 
conformally flat metrics the Dirac equation can be obtained 
from a usual QW with a spacetime dependent coin 11241 1 — 
besides breaking away from general covariance. Yet, ioiiiii 
was no doubt an inspiration : their model can be recovered a 
specific instance of the Paired QW described in this paper, by 
a spacetime grouping as described in Fig. 

Another approach, recently pursued by Succi et al. E2l . is 
within the framework of lattice discretization of the relativis¬ 
tic quantum wave equation (quantum lattice boltzmann ll25ll l. 
The key observation is that the mass term can be recovered 
by extending the neighbourhood of the dynamical map. How¬ 
ever, the existence of a parametrization for the unitary evolu¬ 
tion in terms of the metric, has remained an open question, 
which is solved by this analysis. 

Surprisingly, the slight overgenerality of Eq. Q with re¬ 
spect to the (1 -F 1) curved Dirac equation, just matches some 
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terms arising as (1 + 1) projections of the (2 + 1) curved Dirac 
equation. This suggests a possible generalization to (2 + 1) 
dimensions, through operator splitting, which is the subject 
of current work. The extension to (3 + 1) dimensions re¬ 
mains another interesting question, since gamma matrices of 
the Dirac equation then become four dimensional. In the ter¬ 
minology of this paper this means having to deal with Paired 
QWs having d = 2, so that the internal degrees of freedom of 
the field is four-dimensional. Most of our techniques should 
carry through, except for the characterization of the family 
of solutions to the constraints. Another challenging problem 
is the study the underlying symmetries of the discrete model. 


e.g. by making explicit some form of discrete general covari¬ 
ance along the same lines as 0. 
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Appendix A: Calculation of the first order expansion of the discrete model 

In this section we prove Eq. (|20ll, which we reproduce here: 
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Recall that we want to expand 
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where 

G = E\t + 2, x)W'{t, x){P' © P)[E{t, x-2)® E{t, x + 2)). (A3) 

The first order expansion of the encoding and of the walk is, by definition, 

E{t, x) = x) + eiE^^\t, x)E(t, x) + O(e^) (A4) 

W'{t, x) = x) + x)W{t, x) + O(e^), (A5) 

hence, to first order in e, the operators in ( |A3| ) expand to 

E^t + 2,x)c:i + e (2dtE^°^'f - iEE^^^'f^ (A6) 

W'(t, x) ~ W (A7) 

E{t, x - 2) © E{t, x + 2)- - 2ed^E^°'> + ieE^°^E^ © + 2ed^E(°'^ + leE^^^E^ (A8) 

(A9) 


where in the right hand side all operators are evaluated at (f, x). Recall that the first order expansions of the output and input are 
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(AlO) 


We shall use the identities 


(P' © P){E © E){v © u) = XEv 
(P' © P)(P © E){-v © u) = XZEv 


(All) 

(A12) 


valid for any v G C^, because P' (resp. P) are the projections onto the primed (resp. non-primed) coordinates; in matrix form. 


/O 0 1 0\ p _ A 0 0 0\ 

i^oooij’ 


Next we plug the previous expansions into (|A2[). Collecting all the terms of first order in s. 
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Next we use the zeroth order condition (cf. @), namely p(o)tVp(o)xp(o) = J © P, so that 
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and we get the desired result. 





































9 


Appendix B: General form of B 


Since B must be hermitian, cf. ( |21a| l, then Bi and B 4 are hermitian. Since it is also unitary, then it must square to the identity. 
This implies that the conditions 


Bl + bIB2 = Id2 

(Bla) 

5| + B 2 BI = Id2 

(Bib) 

and 


B 2 B 1 -h B^B 2 — 0 

(B2a) 

B 1 B 2 “h B^B^ = 0 

(B2b) 

must hold. Note also that B must have a complete set of orthonormal eigenvectors, eigenvalues ±1, and it shall be traceless. 


because it is is similar to Z. 

First, we parametrize the block B 2 . Consider the spectral decomposition of Bi = VDV\ D = diagjdi, ^2}. From the 
first of conditions ( |Blal l, we have that di,d 2 G [—1,1], because the square root of the components of Id —D^ are precisely the 
singular values of B 2 , which should be non-negative. Next, we shall find B 2 such that constraint ( |34a| ) is satisfied. The same 
equation also determines U. 

We look for B 2 G such that conditions ( |34a] i and ( |Bla] i are satisfied, namely that 

1. Id+252 is unitary, 

2. B^B2 = Id-Bf. 

To prove our lemma we will use a shortcut provided by the following characterization of matrices with positive definite 
real part. Recall that if A £ C"^", its real and imaginary parts are 3?+ := ^(A + and 5A = respectively. 

Theorem 1 (see ll^Tl l. Let A G C"^". Then, ^A is positive definite if and only if 

/l + icki 


A = T 




(B3) 


V 


1 + ickn 


for some non-singular T and ai,..., a„ £ K. 

Note that, from condition 1 above, 

(Id+25|)(Id+252) =ld^ B 2 + bI + 2(5|52) = 0, 

hence condition 1 is equivalent to 3?52 = —B\B 2 . Recall that A^ A is positive definite for any A £ 
can be applied to — ^2. 

We recall the following parametrization of the U(2) group, namely that 


(B4) 

hence Theorem 


C/(2) = e- 




:0£ [0,27r), a,/?£C, jap + |/3|2 = 1 


(B5) 


Lemma 1. Let Bi G be hermitian, with spectral decomposition Bi = ViDiVi, and eigenvalues di, d 2 G [—1,1] £ 

Assume that B 2 G satisfies the conditions 

1. 3?52 = -BIB2, 


2. bIB 2 = Id-Bl 

Let Xi = Vl — df, and rji G {flt^Vi}’ V 2 G with ijf G [0,7r/2], r/~ G [—7r/2,0], and such that sin rjf = zL\di\for 

any i G {1,2}. Then, 

1. If d\ ^2 (non-degenerate case) then 


B 2 = -Fi 


0 

0 A2e‘''" 




(B6) 
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2. If df = d 2 and rji = 772 , then B 2 = Id. 

3. If d\ = ^2 and rji = — 772 . then 

B 2 = 


(B7) 


for any K G U(2). We remark that any two K 12 such that Ki = K 2 K' for some K' of the form K' = 
cos(0) Id+isin( 0 )(T 2 , will give the same i? 2 - 

Proof From Theoremj^ and condition 1, B 2 can be written as B 2 = —T 2 D 2 T 2 for some non-singular T 2 and D 2 = diag{l -1- 
iai, 1 -I- ia 2 }, cti, 02 € K. Substitution into condition 1 gives 


nT 2 = 


(l-baf) ^ 0 

0 (1 + ai)-' 


(B8) 


Now, let the SVD of T 2 = . Then T 2 T 2 = and using ( |B 8 | ) it is easy to see, using the canonical decomposition 

of unitary matrices (cf. ( |B5| l), that we must have one of the following cases: 

1. If af f a|, then for some 6 * 1,02 € [0, 27r), either 


(a) S = 


_ I yi+o 


0 


(b) s = 


_ ( ^\+a\ 


0 


0 

1 

■y/l+c 


and V = 


and V = 


0 

,je2 0 


0 

0 


Hence, B 2 =-W 


/ l+ia2 
/ T+Tf 


Hence, B 2 =-W i 


l+iai 


l + ia:2 
T+aF, 






Note that in either case, we can write B 2 = —H^diag W\ where cr : {1, 2} —>■ {I, 2} is a 

permutation. It is easy to check that condition 1 is indeed satisfied. Next, substitution into condition 2 gives 


1 — d\ 

0 


0 

1-di 


= K 


l+“o 

0 


(1) 


k\ 


(B9) 


1+a^ 


( 2 ) . 


with K ■= V^W € U(2). Again, we use the canonical form, cf. (HD, to find that K is diagonal or antidiagonal, 

(j2 

with two independent phases. Introducing back W into B 2 , in either case we have i = 1,2, hence 

I + ia* 


Oti = ± SO 


1 + af 


s/l- d} 1 - df ± i|di|^ = \J\- d}e"^^, 


(BIO) 


provided that cos 77 ^ = s/l — dif, and sin 77 ^ = ±|(ii|. This proves part 1. 


2. If af = a|, then E = 


, ^ Id 2 , and V € U(2) is arbitrary. Thus, we can write B 2 = — ,} 2 K 


1 + iai 0 
0 1 ± iai 




for some K G 17(2). Substitution into condition 2, gives af = df/(l — df), and proceeding as in ( |B10| l, we obtain the 
claim. 


□ 

Next we characterize B 4 . We assume that the relevant constraints fron ( |Bla| i-( |B2b| ) are satisfied, namely Bf + b\B 2 = Id 2 
and that i3| + B 2 BI = Id 2 . 

Lemma 2 . In the hypothesis of above, 

1. If d\ f 02 (non-degenerate case), then B 4 = —Bi. 

2. If di = 02, then B 4 = —Bi. 
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5. Ifdi = —d 2 , then 


Bi = diKa^K\ (Bll) 

for any K € U(2). We remark that any two K 12 € U{2) such that Ki = K 2 K' for some K' of the form K' = 
cos(0) Id +i sin(0)(T2, will give the same B^. 

Proof From ( |34a| l, B 2 is normal, then from ( |Bla| l and ( |Blb| l we have that B\ = B\. Since ^4 is hermitian consider its spectral 
decomposition, i ?4 = WD 4 W\ W G U{2). Then, Dj = KDfKf where K := Using the canonical form ( |B5| l, we 

hnd that 

1. If dl 7 ^ ^ 2 ’ '^hen K is either diagonal or anti-diagonal, with arbitrary phases. In either case we obtain i ?4 = 
Uidiag{±di, ±d 2 }U]^, but since we must have Tri 34 = —Tri?i, we shall take —di, —d 2 - Hence, i 34 = —Bi. 

2. If d\ = df, then K G U{2) is arbitrary, and we have B 4 = ViiT^diag{±di, ±di\KVi. If di — d 2 , then Tri 34 = —2di, 
so B 4 = —di Id 2 = —Bi. If di = —d 2 , then TrB 4 = 0, and we can take di, —di or —di, di. 


□ 







